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Abstract. - The phase diagram of binary mixtures of particles interacting via a pair potential 
of parallel dipoles is computed at zero temperature as a function of composition and the ratio 
of their magnetic susceptibilities. Using lattice sums, a rich variety of different stable crystalline 
structures is identified including AmBn structures. [A (B) particles correspond to large (small) 
dipolar moments.] Their elementary cells consist of triangular, square, rectangular or rhombic 
lattices of the A particles with a basis comprising various structures of A and B particles. For 
small (dipolar) asymmetry there are intermediate AB2 and A2B crystals besides the pure A and 
B triangular crystals. These structures are detectable in experiments on granular and colloidal 
matter. 
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While the freezing transition and the correspond- 
ing crystal lattice in one-component systems is well- 
Understood by now [1,2], binary mixtures of two different 
particle species exhibits a much richer possibility of differ- 
ent solid phases For example, while a one-component hard 
sphere system freezes into the close-packed face-centered- 
dubic lattice [3], binary hard sphere mixtures exhibit a 
huge variety of close-packed structures depending on their 
diameter ratio. These structures include ABn superlat- 
tices, where A are the large and B the small spheres, with 
n = 1, 2, 5, 6, 13. These structures were found in theoreti- 
cal calculations [4], computer simulations [5,6] and in real- 
space experiments on sterically stabilized colloidal suspen- 
sions [7,8]. Much less is known for soft repulsive interpar- 
ticle interactions; most recent studies on crystallization 
include attractions and consider Lennard- Jones mixtures 
[9,10] or oppositely charged colloidal particles [11-13]. 

In this letter we explore the phase diagram of a binary 
mixture interacting via a soft repulsive pair potential pro- 
portional to the inverse cube of the particle separation. 
Using lattice sums, we obtain the zero-temperature phase 
diagram as a function of composition and asymmetry, i.e. 
the ratio of the corresponding prefactors in the particle- 
particle interaction. Our motivation to do so is threefold: 

i) First there is an urgent need to understand the ef- 
fect of softness in general and in particular in two spatial 



dimensions. The case of hard interactions in two spatial 
dimensions, namely binary hard disks, has been obtained 
by Likos and Henley [14] for a large range of diameter ra- 
tios. A complex phase behavior is encountered and it is 
unknown how the phase behavior is affected and controlled 
by soft interactions. 

ii) The model of dipolar particles considered in this let- 
ter is realized in quite different fields of physics. Dipo- 
lar colloidal particles can be realized by imposing a mag- 
netic field [15]. In particular, our model is realized by 
micron-sized superparamagnetic colloidal particles which 
are confined to a planar water-air interface and exposed 
to an external magnetic field parallel to the surface nor- 
mal [15-19]. The magnetic field induces a magnetic dipole 
moment on the particles whose magnitude is governed by 
the magnetic susceptibility. Hence their interaction poten- 
tial scales like that between two parallel dipoles with the 
inverse cube of the particle distance. Binary mixtures of 
colloidal particles with different susceptibilities have been 
studied for colloidal dynamics [20], fluid clustering [21,22], 
and the glass transition [23]. A complementary way to 
obtain dipolar colloidal particles is a fast alternating elec- 
tric field which generates effective dipole moments in the 
colloidal particles [24]. This set-up has been applied for 
two-dimensional binary mixtures in Ref . [25] . In granular 
matter, the model has been realized by mixing millimeter- 
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sized steel and brass spheres [26] which are placed on a hor- 
izontal plate and exposed to a vertical magnetic field such 
that the repulsive dipolc-dipolc interaction is the lead- 
ing term. Stable triangular AB2 crystalline lattice were 
found [26] . Layers of dusty plasmas involve particles whose 
interactions can be dominated by that of dipolcs [27-29]. 
Other situations where two-dimensional mixtures of par- 
allel dipolcs are relevant concern amphiphiles (confined 
to a monomolccular film at an air- water interface [30]), 
binary monolayers [31-33], /erro/?wic! monolayers [34] ex- 
posed to a perpendicular magnetic field, or thin films of 
mo/ecu/ar mixtures (e.g. of boron nitride and hydrocarbon 
molecules) with a large permanent dipolc moment [35]. 
Hence, in principle, our results can be directly compared 
to various experiments of (classical) dipolar particles with 
quite different size in quite different set-ups. 

iii) It is important to understand the different crys- 
talline sub-structures in detail, since a control of the col- 
loidal composite lattices may lead to new optical band-gap 
materials (so-called photonic crystals) [36] to molecular- 
sieves [37] and to micro- and nano-filters with desired 
porosity [38]. Nano-sieves and filters can be constructed 
on a colloidal monolayer confined at interfaces [38] . Their 
porosity is directly coupled to their crystalline structure. 
For these applications, it is mandatory to understand the 
different stable lattice types which occur in binary mix- 
tures. 

As a result, we find a variety of different stable com- 
posite lattices. They include AmBn structures with, for 
instance, n= 1,2,4,6 for m = 1. Their elementary cells 
consist of (equilateral) triangular, square, rectangular and 
rhombic lattices of the A particles. These are highly deco- 
rated by a basis involving either B particles alone or both 
B and A particles. The topology of the resulting phase 
diagram differs qualitatively from that of hard disk mix- 
tures [14]. For small (dipolar) asymmetries, for instance, 
we find intermediate AB2 and A2B structures besides the 
pure triangular A and B lattices which are absent for 
hard disks. Our calculations admit more candidate phases 
than considered in earlier investigations [39] where two- 
dimensional quasicrystals were shown to be metastable. 
We further comment that we expect that colloidal glasses 
in binary mixtures of magnetic colloids [23] are metastable 
as well but need an enormous time to phase separate into 
their stable crystalline counterparts. 

The model systems used in our study are binary mix- 
tures of dipolar particles made up of two species denoted as 
A and B. Each component A and B is characterized by its 
dipole moment m.A and ms, respectively. The particles 
are confined to a two-dimensional plane and the dipole 
moments are fixed in the direction perpendicular to the 
plane. Thereby the dipole-dipole interaction is repulsive. 
Introducing the ratio m = tub I m a of dipole strengths 
rriA and ms, the pair interaction potentials between two 
A dipolcs, a A- and B-dipole, and two B-dipoles at dis- 



tance r are 

VAA{r)= Vbv3(r), VAB{r) = Vbm(p(r), 

V^si3(r) = yomVW, (1) 

respectively. The dimensionlcss function ^p(r) is equal 
i'^/r^, where £ stands for a unit length. The amplitude 
Vo sets the energy scale. 

Our task is to find the stable crystalline structures 
adopted by the system at zero temperature. We consider 
a parallelogram as a primitive cell which contains ua A- 
particles and riB i?-particles. This cell can be described 
geometrically by the two lattice vectors a = a (1,0) and 
b = a7(cos6',sin6'), where 9 is the angle between a and 
b and 7 is the aspect ratio (7 = |b|/|a|). The position 
of a particle i (of species A) and that of a particle j (of 
species B) in the parallelogram is specified by the vectors 
= {xf-^yf^) and = {x^,y^), respectively. The total 
internal energy (per primitive cell) U has the form 

J=A,Bi,j=l R 

+ T.Y.Y.^^B{\vt^rf+K\), (2) 

i = l j = l R 

where R = ka + Ih with k and I being integers. The 
sums over R in Eq. 2 run over all lattice cells where the 
prime indicates that for R = the terms with i = j are 
to be omitted. In order to handle efficiently the long- 
range nature of the dipole-dipole interaction, wc employed 
a Lekner-summation [40,41]. 

We choose to work at prescribed pressure p and zero 
temperature (T = 0). Hence, the corresponding ther- 
modynamic potential is the Gibbs free energy G. Ad- 
ditionally, we consider interacting dipoles at composition 
X := TiB/inA + nB), so that the (intensive) Gibbs free en- 
ergy g per particle reads: g = g{p, m, X) = G/ {tia + tlb)- 
At T = 0, (7 is related to the internal energy per particle 
u — U/ (jiA + riB) through g — u+p/ p, where the pressure 
p is given hy p = p'^{du/dp), and p — {ua + nB)/|a x b| 
is the total particle density. The Gibbs free energy per 
particle g has been minimized with respect to 7, 9 and 
the position of particles of species A and B within the 
primitive cell. To reduce the complexity of the energy 
landscape, we have limited the number of variables and 
considered the following candidates for our binary mix- 
tures: AiB, A3B, A2B, A4B2, A3B2, AB, A2B2, A2B3, 
AB2, A2Bi, AB3, AB4 and ABq. For the ABq case we 
considered a triangular lattice formed by the A particles. 

The final phase diagram in the (m, X)-plane has been 
obtained by using the common tangent construction. The 
dipole-strength ratio m can vary between zero and unity. 
A low value of m (i.e., close to zero) corresponds to a 
large dipole-strength asymmetry, whereas a high one (i.e., 
close to unity) indicates a weak dipole-strength asymme- 
try. Our calculations show that all the mixtures, except 
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Fig. 1: The stable binary crystal structures and their primitive cells. The red dark (green light) discs correspond to A (S) 
particles. 
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Table 1: The stable phases with their Bravais lattice and their 



basis. 
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AB^i and A4B2, are stable. Their corresponding crys- 
talline lattices are depicted in figure 1 and the nomen- 
clature is explained in Table 1. For the one component 
case [X — (pure A) and X — 1 (pure B), see figure 2], 
we found an equilateral triangular lattice T(A) and T(i3), 
respectively, as expected (see figure 1). 

The most relevant and striking findings certainly con- 
cern the phase behavior at weak dipole-strength asymme- 
try (0.5 < TO < 1), see figure 2. Thereby, the only stable 
mixture ^1-62 over such a large range of to corresponds 
to the ("globally" triangular) phase T{AB2) (see figure 1 
and figure 2). This is in strong contrast to what occurs 
with hard disk potentials [14], where no mixture sets in 
at low size asymmetry. At sufficiently low dipole-strength 
asymmetry (m > 0.88), see figure 2, the mixture A2B, 
that also corresponds to a globally triangular crystalline 
structure [namely T{A2B), see figure 1], is equally stable. 
The stability in the limit m ^ 1 of those globally tri- 
angular structures are fully consistent with the fact that 
one-component dipolar systems are triangular. 

In the regime of strong dipole-strength asymmetry 
(0.06 < m < 0.5); sec figure 2, the stability of the com- 
position X = 1/2, corresponding to the mixtures AB and 
A2B2, is dominant and the phase diagram gets richer in- 
volving all the different structures [except T{AB2)] shown 
in figure 1. More specifically, for the composition X = 1/2, 
we have two phases S{AB) and Ilh{A)AB2. The transi- 
tion between these two phases is continuous as marked by 
a symbol # in figure 2. For X = 2/3, many stable phases 
emerge as depicted in figure 2. In the i?-rich region at 
large asymmetry, the stability will involve many different 



structures which are probably not considered here. There- 
fore we leave this region open, see the gray box in figure 2. 
Below X = 2/3, at large asymmetry (to < 0.2), the true 
phase diagram will also involve a very dense spectrum of 
stable compositions, as suggested by the already many sta- 
ble compositions (see figure 2), which are not among the 
candidate structures considered here. This feature is very 
similar to the behavior reported in hard disk mixtures [14], 
where a continuous spectrum of stable mixtures is found 
for X < 2/3 at high size asymmetry. In the limit to —^ 0, 
a triangular lattice for the A particles will be stable with 
an increasingly complex substructure of B particles. 

In conclusion, the ground-state phase diagram of a 
monolayer of two-dimensional dipolar particles shows a 
variety of different stable solid lattices. The topology of 
the phase diagram is different from that of hard disks. 
Whereas short-ranged interactions lead to a phase separa- 
tion into pure A and B crystals at low asymmetries, there 
are two intermediate A2B and AB2 mixtures for softer 
interactions. This explains the experimental findings of 
Hay and coworkers [26] who found an AB2 crystal struc- 
ture in millimeter-sized steel and brass spheres [26] which 
does not occur to be stable for hard particles. A further 
more quantitative experimental confirmation of our theo- 
retical predictions are conceivable either in suspension of 
magnetic colloids or for binary charged colloidal suspen- 
sions [42] confined between two parallel glass plates [43] 
or for any other situation where two-dimensional dipolar 
particles are involved. 

We finish with a couple of remarks: First, based on 
the present studies it would be interesting to study the 
behavior of tilted dipoles where anisotropics and attrac- 
tion play a significant role [44]. Our data may also serve 
as a benchmark to perform further studies on melting of 
the composite crystals and crystal nucleation out of the 
melt in two spatial dimensions. The extension to one- 
component bilayers [45, 46] made up of dipolar particles 
would certainly be relevant. It would also be interesting 
to apply the method of evolutionary algorithms [47] to the 
present problem in order to increase the basket of candi- 
date phases. 

We thank Christos Likos, Hans- Joachim Schope and 
Thomas Palberg for helpful discussions. This work has 
been supported by the DFG within the SFB-TR6, Project 
Section Dl. 
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